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The gravity-scalar field system in spherical symmetry provides a natural setting for exploring 
gravitational collapse and its aftermath in quantum gravity. In a canonical approach, we give 
constructions of the Hamiltonian operator, and of semiclassical states peaked on constraint free 
data. Such states provide explicit examples of physical states. We also show that matter-gravity 
entanglement is an inherent feature of physical states, whether or not there is a black hole. 

I. INTRODUCTION 

Black hole thermodynamics, cosmology and more generally unification, have all provided strong impetus for de- 
veloping a theory of quantum gravity. There appears to be agreement on at least two features such a theory should 
have regardless of the details of the approach. These are background (or metric) independence and fundamental 
discreteness. It may be that other ideas such as holography arise as a consequence of the first two but this is far from 
clear at the present level of understanding. 

Progress in quantum gravity has often been limited to models where a complete quantization can be performed, such 
as the work on mini-superspace cosmology in the early seventies [1], and its later incarnations in the Hartle-Hawking 
approach [2] and loop quantum cosmology [3]. These are all quantum mechanical systems arrived at by the assumption 
of spatial homogeneity of spacetime. The extent to which such models can reveal insights into quantum gravity is not 
clear [4] although recent developments on singularity avoidance may be one such feature. More generally the hope is 
that such smaller systems will ultimately become absorbed in a full theory of quantum gravity in the same way that 
the Bohr atom has in quantum electrodynamics. 

Going beyond mini-s\ipcrspace models significantly enlarges the problem to a true field theory of a constrained 
system. The first step in this direction is to a two-dimensional field theory, where the metric and matter variables 
depend on time and one spatial coordinate. There are a few models of interest in this category, namely the Gowdy 
cosmology, cylindrical gravitational waves and the asymptotically flat gravity-scalar field theory. The quantum theory 
of the first two models have been studied in the simplified approximation of one local degree of freedom [5, 6]. These 
are vacuum models so they are less interesting physically than the last one, which has the potential to reveal much 
about gravitational collapse and black hole formation in quantum gravity [7-9]. Furthermore this model is the natural 
next step beyond just the quantum mechanics of the Schwarzschild black hole, on which much has been written [10-15] . 

This system has a true Hamiltonian that forms a part of its asymptotic Poincare symmetries. The model provides 
the setting for two important physical scenarios: numerical studies of gravitational collapse [17, 18] and the so-called 
information loss paradox, whose origin involves an assumption on the large scale structure of an evaporating black hole 
based on semiclassical physics [16]. It is concerned with a collapsing low-entropy (or even pure-state) matter that forms 
a black hole, which subsequently evaporates within a finite time, resulting in a highly-cntropic Hawking radiation 
[19]. If the correlations between the inside and outside of the black hole are not restored during the evaporation 
process, then the corresponding increase in entropy is interpreted as a lost "information" . Various arguments have 
been advanced to justify or dismiss this phenomenon [20], as well as possible consequences of either option [21]. 
There is the possibility that the existence of a quantum theory of gravity automatically removes the information 
loss "paradox", since a complete quantization of the true Hamiltonian would result in the unitary evolution of the 
combined matter-gravity system [22] , while the entropy flow between subsystems can be interpreted with the help of 
quantum information theory [23]. 

In this paper we build on recent developments [7, 8] aimed at obtaining a compete quantum theory of the asymptot- 
ically flat gravity-scalar field theory in spherical symmetry. This earlier work used a gauge appropriate for Painleve- 
GuUstrand coordinates (because of their regularity at the horizon), and provided a construction of the Hamiltonian 
and null expansion operators. Not addressed in these works are the problems of finding physical states and addressing 
their dynamics. This is what we study in this paper. 

The outline of this paper is as follows. We first review the model and its quantization. We then give a prescription 
for constructing semiclassical physical states that are peaked on classical data, and give a new form of the Hamiltonian 
operator. (An earlier version of this operator utilized a Dirac-like trick to write a square root operator.) Using these we 
show that physical states exhibit matter- gravity entanglement. We argue that this may be a robust model independent 
feature of non-locality in quantum gravity 
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II. GRAVITY-SCALAR FIELD MODEL 

Our starting point is the Arnowitt-Deser-Misner (ADM) hamiltonian formulation for general relativity. The phase 

space of the model is defined by prescribing a form of the gravitational phase space variables Qab and 7r°^, together 
with fall-off conditions for these variables, and for the lapse and shift functions N and N'^. The ADM 3+1 action for 
general relativity minimally coupled to a massless scalar field is 



S = J d'xdt [TT^^-gab + P^^ -NH- N'^Ca 

The pair P^) are the scalar field canonical variables, and 



(1) 



+87rG' (y^Pl + Vq<l"''da<pdbcp^ ~ (2) 
Ca = i?c< - 87rGP05„0 ~ 0, (3) 

where tt = Tr°-^qab- 

This action (together with the boundary terms, see e.g. [8, 24]) is well-defined and determines the fall-off conditions 
on canonical variables. Below (Sec IIA) we outline the structure that results from imposition of spherical symmetry 
and use of the fiat-slice partial gauge- fixing. (This is a summary of the work in [8].) The reason for utilizing this 
gauge is that the horizon is not located at a coordinate singularity. In Sec IIB we present a family of initial data that 
will illustrate, in a later section, a construction of semiclassical physical states of the model. 

A. Hamiltonian theory in the flat slice formalism 

The reduction to spherical symmetry utilizes an auxiliary flat Euclidean metric Cat and unit radial normal s" = x^/r, 
where = Cabx'^x^- The parametrization of the reduced phase space we use is given by the matrices 

qab = A(r,t)^SaSfcH ^ — {eab-^aSb), (4) 

P4r,t) .,^r^PR{r,t) 
2A(r,t) 4R{r,t) 



= s'^s' + rJ:" ' (e"^" - s'^s"), (5) 



As a line element the spatial metric is therefore 

ds'^ = A^{r,t)dr^ + R{r,tfdn^, (6) 

where the solid angle d^l^ arises from the second term in these expressions. 

The Painleve-GuUstrand (PG) coordinates are those where equal coordinate time slices are spatially flat [8, 24]. 
However it is sufficient to use the partial gauge flxing A = 1 to obtain the feature of non-singular coordinates at the 
horizon, which is the feature of PG coordinates we desire. This is what is summarized in the remaining part of this 
section. 

Substituting (5) into the ADM 3 + 1 action with minimally coupled scalar fleld leads to the reduced action 

= ^ j dtdr [PrR + FaA + P^<^) 
~4^/ dtdriNH + N'-Cr) 

dt{N-APA)\r=.o, (7) 



where N,N'^ are the lapse and radial shift functions, and the reduced Hamiltonian and (radial) diffeomorphism 
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constraints H and C" are 



H = 



i?2A 



\ {P^Kf-\{P^K){PnR) 



+ -r^ \2RR!'K - 2RR'A' - + AR'- 



+ 



i?2 



2Ai?2 2A 
Pi^i?' - AP; + 



.0, 
0. 



(8) 
(9) 



These constraints are first class with an algebra that is similar to that for the full theory. 

We note that since this reduction to spherically symmetry is unusual in that the usual ADM mass integral vanishes 
due to the flat slicing condition. It is therefore important to ensure that the reduced action in functionally diffcrcntiablc 
with a consistent set of fall-off conditions on the phase space variables, and on the lapse and shift functions. This 
analysis has been carefully done by one of the authors in Ref. [8], where these details are explicitly given. The 
asymptotic conditions on the phase space variables are 



R 

Pr 
A 

Pa 



and those on the lapse and shift functions are 

N 



Ar-^/y2 + Oir-^-'), 
l + 0(r-3/2-^), 



= Ar 



-1/2 



0{r 



-l/2-eN 



(10) 

(11) 
(12) 

(13) 



where A,B,C are constants. The leading order terms are solutions of the constraints, so there is no logarithmic 
divergence of the action at this order. The mass formula is the surface term in Eqn. (7); the mass parameter arises in 
this integral through the momentum Pa- It is readily verified that if the scalar field is set to zero the Schwarzschild 
solution results [8]. 



B. Partial gauge fixing 

We next impose the gauge choice A = 1, which corresponds to a step toward flat slice coordinates. With this 
condition imposed, the Hamiltonian constraint is solved (strongly) for the conjugate momentum Pa as a function of 
the phase space variables. This gives 

Pa = PrR + ^{PrRY - X, (14) 

where 

a: = 16P2(2PP"- 1 + P'2) + I6p2if_^^ (15) 

and 

This effectively eliminates the conjugate pair (A, Pa) in favour of the remaining phase space variables. We must 

however obtain the consequence of this for the lapse and shift functions in the standard way: the evolution equation 
for A [8] and the requirement that the gauge A = 1 be preserved under evolution leads to the equation 

N = , ^ ' . (17) 

V(PflP)2 - X 
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We note that this relation docs not lead to any restriction on the phase space variables, since it is merely a relation 
between the lapse N and and the shift TV, neither of which is absolutely fixed at this stage. We also note that on the 
classical space of solutions, the argument of the square root is positive; in quantum theory, it is in principle possible 
that fluctuations may make this negative. This however is an issue that is present in general in the ADM variable 
approach to quantum gravity because the Hamiltonian constraint contains a factor of y/q, which may fluctuate such 
that the argument is negative. The issue hinges on how the corresponding operator is defined. We discuss this below 
where we set up the Dirac quantization problem. 

The gauge A = 1 reduced the radial diffeomorphism constraint to 

Cr = -P; + PrR' + P^(j)' ~ 0, (18) 

with Pa given by (14) above. We note that using this constraint the square root in the latter equation can be written 
as 

^{PrR^ -X= j {PrR' + P4,(t>') - PrR (19) 
Jo 

To summarize this section, the partially gauged fixed theory is prescribed by the phase space variables {(j), P^) and 
(P, Pr) and the reduced Hamiltonian 

/■CX) 

= / [{N-^yPK + N-'iPRR! + P4,<l>')]dr 

JQ 

/•OO 

+ N^iPnR' + P^(t>') dr, (20) 
Jo 

where the surface term in the reduced action (7) has been written as bulk term and combined with the remaining 
radial diffeomorphism constraint. Together with the asymptotic conditions given above, the variational principle is 
well defined. This is the reduced system we study in the rest of the paper. 

Finally we note that we have not yet fully fixed the gauge to the flat slice case since we have not imposed the 
condition R = r. We prefer to retain P as a dynamical variable for developing the quantum theory, since full gauge 
fixing leads to a reduced Hamiltonian for the scalar field degrees of freedom that is non-local [8], and so poses a 
problem for quantization. 



C. Classical data 



Finding classical initial data sets is necessary both in numerical calculations and for the construction of semi-classical 
physical states. As one of our aims in this paper is to present a concrete realization of the latter, it is useful to see 
what form such data takes before we proceed to a construction of such states. The data should be an asymptotically 
fiat solution of the remaining constraints prescribed by functions (^(r),P^(r) and P(r) and PR{r). The family of 
data that we exhibit below depends on two parameters. It is regular at the coordinate origin and may have trapped 
surfaces, depending on the parameter values. 

We shall see, after describing the quantum theory, that semiclassical states corresponding to data with trapped 
surfaces such as this, describe a quantum black hole state. This occurs because states peaked on classical constraint 
free data satisfy the condition that the expectation value of the constraint operator vanishes to leading order in a well 
defined expansion (see Section V below.) 

With the partial gauge fixing A = 1 already imposed as described in Sec. H.B, we fix gauge fully by imposing 
P = r. This simplifies the constraint to 

Cr = -P^ + PR + P^(l>' = 0, (21) 

with 

Pa = Pijr + V(Pflr)2 - X, (22) 

where X now reduces to 

X = 16P^(2PP" - 1 + P'2) + 16R^H^ = l6r^H^. (23) 
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To find some explicit solutions of interest, we make the ansatz (p = 0. The constraint then becomes a relation between 
Pr and P^, so one of these may be chosen freely. 
To obtain an explicit class of solutions, let us set 

16H^ = h\r)Pl{r), (24) 

for some function h. Then (22) becomes 

Pa = Pflr(l + Vl - h?) = PRr{l + g), (25) 
and the auxiliary function h = \J\ — satisfies 

< /i^ < 1, lim = 0. (26) 

r— )-oo 

The constraint is now 

P'j,r + [PRvg)' = 0, (27) 

which can be rewritten as 

{Pnr)' -Pr + {PRvg)' = 0. (28) 
Finally, this may be put into an integrated form by setting 

Pr = n'fl (29) 

to give 

Thus, given suitable functions g, wc can find Hr and hence Pr, followed by P^ from Eqn. (24). 
We note that the energy density can be written as 



,iirPr n 



16//^ = 2^-^. (31) 

The asymptotic fall-off conditions (above) require Pr ~ l/v^, which implies ~ 2PRr and <? — >■ 1. If one wishes to 
have a regular data at r = 0, then the expansion IIr ~ r" leads to 

< - - 1 < 1, (32) 
a 

which bounds the power to 

1/2 < a < 1, (33) 

while the regularity of imposes a = 1, so for r we get IIr ~ cr 

To obtain an explicit asymptotically flat solution let us consider the choice 

l+ff=(l-iexp(-a/a;2))-i, (34) 

which results in 



Ur = cr exp 

where the integral exponential function is deflned as 



(35) 



/•oo t 

Ei{z) = J —dt. (36) 
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FIG. 1: Expansion of the outgoing null geodesic with (solid line) and without (dashed line) trapped surfaces on the initial 
hypersurface. The solid line corresponds to o = 1, C = 2 and the dashed line to o = 1, C = 6). 



The asymptotic behavior of JIr is as desired, 

\Ir ~ cyae''/^y/r, r ^ oo, 

where Euler's ganiina is 7 = 0.57 . . ., and 

Hji cr, r ^ 

From the asymptotic fall-off condition 

PR = ylr-V2/2 + 0(r-i-^), 

where A is related to the ADM mass as A = 4a/2M [8], we see that cv^c"'/"' = A. 

The next step in the investigation of the classical problem is to check for the presence of apparent horizons. In the 
gauge we adopted the expansions of outgoing and ingoing null geodesies are given as [7] 



(37) 
(38) 
(39) 



that reduces for the initial data surface to 



e± = ±4:RR' - Pa, 



e± = ±4r - Hr 



(40) 



(41) 



By varying parameters of Ur it is possible to generate the regular initial data with or without trapped surfaces, as 
shown on Fig. 1. 



III. QUANTIZATION 

In this section we describe an approach to constructing the quantum theory of the system described above. We 
will be working with the partially gauge fixed theory described in Sec. II. B, where we have the pair of phase space 
variables (i?, P^) and (</), P^), and the remaining radial diffcomorphism constraint. The quantization we use is a 
type of polymer quantization where field configuration and translation operators are quantized, but there is no direct 
quantization of field momenta. Using the basic operators, we shall see that there are well defined prescriptions for 
writing the constraint and the Hamiltonian operators. 



A. Hilbert space and basic operators 



To describe the polymer quantization [26, 27] (which may be viewed as the "dual" to that used in loop quantum 
gravity), we begin with the basic variables 

poo 

Rf= drf{r)R{r), [/^(Pij) = exp (ieP^) , (42) 
Jo 
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where /(r) is a smearing function and e is a real dimensional constant. These satisfy the canonical Poisson bracket 

{Rf, U,{PH{r))} = ief{r)Ux{PRir)). (43) 

(In defining these observable we have at our disposal a constant density that ensures the correct density weights in 
the integral for Rf and in the exponent of 11^. This density is the phase space variable A which was set to unity when 
we performed the partial gauge fixing A = 1. ) 

We use similar definitions for the variables made from (f) and P^, ie. 

/•oo 

<^/= / dr/(r)0(r), ^^(P^) = exp (mP^) . 

Jo 

The Poisson bracket of these matter variables mirrors that of the metric variables in (43) . 

The parameters e and k have the physical dimensions necessary to make the exponents in the corresponding U 
variables dimensionless. This scale will carry over to quantization. Since we are dealing with a quantization of gravity 
and matter, it is natm-al to suppose that that both parameters are related to the Planck scale. In the following we 
will keep them separate. 

The Poisson bracket of the metric variables is realized as an operator relation on a Hilbert space spanned by the 
basis states 

|ai,a2,--- ,a„) (44) 

where the real numbers Oj represent values of the configm-ation variable R at the radial points fj. The points 
provide a set where the phase space variables are sampled. The selection of space points rj may be thought of as a 
lattice (or graph in the language used in loop quantum gravity). Since there is a basis state of this type for every 
countable set of points, the Hilbert space is non-separable. In the following it will be convenient to work with a fixed 
and uniformly spaced set of points, although this is not necessary. 
The inner product is 

{a[,a'2, ■ ■ ■ ,ajai,a2, ■■■ ,a„) = Sa'^^i ' ' ' Sa'^,a„ (45) 

if two states are associated with the same lattice points; if not the inner product is zero. This inner product is 
background independent in the same way as for example the inner product for the Ising model; the difference is that 

for the latter there is a finite dimensional space of spins at each lattice point. 

The configuration and translation operators are defined by the following expressions: 

Rf\ai,a2, - ■ ■ cin) ■=^aif{ri)\ai,a2,---an), (46) 



Ue{PR{rk))\ai,a2,---an) ■= |ai,a2,--- , afe - e, • • • a„). (47) 

It is readily verified that the commutator of these operators as defined provides a faithful realization of the corre- 
sponding Poisson bracket. Since the U operators are realized here as ladder operators for the field excitations, the 
parameter e represents the discreteness scale in field configuration space. 

This representation is one in which the momentum operator does not exist. There is however an alternative e 
dependent definition of momentum using the translation operators (47), given by 

Phirk) := 2^ {UeiPnirk)) - UKPR{rk))) , (48) 

which will be used in the definition of the Hamiltonian operator. 

The representation for the matter variables is similar to that defined above for the metric variables. We write the 
basis states of the matters sector as 

\bi,b2,--- ,bn) 

With the inclusion of matter the kinematical Hilbert space is the tensor product of geometry and matter Hilbert 
spaces, with basis 

|ai,...,ajv;6i,...,6]v), (49) 

^ ^ ^ -y ^ 

gravity matter 
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With the quantization of the basic variables in hand, we now turn to defining the composite operators necessary to 
write the constraints and the Hamiltonian. To this we first need to define a locahzed field operator from Rf and (j)f. 
A localized field may be defined by taking for example /(r) in eqn. (42) to be a Gaussian 

G(r,rfe,(T) = e-^'--'-")'/'^', (50) 

(or a smooth function of bounded support) which is sharply peaked at a radial point r^- For illustration we will work 
with the Gaussian with the understanding that its width a is such that the function is effectively zero at all lattice 
points except where it is peaked. This is easiest to see with a uniform lattice such that cr <C 1 in Planck units. With 
this in mind we write 

^G{rk) = ^k, (51) 

with a similar expressions for the scalar field </>. We also have the following action of the basic operators (where we 
have included explicitly the Planck length): 

Rk\ai, . . . , Oat; 6i, . . . , 6jv) = 2l%ak\ai, . . . , oat; 6i, . . . , 6jv), (52) 



(f>k\ai, ■ ■ ■ ,aN] h, - ■ ■ ,^jv) = 2Zp6fe|oi, . . . ,aN;bi, . . .,bN)- 
The field translations operators f/fe(e) = e'*^K'= and Vk{K) = e*'*^*'= act as 

Uk{^)\ai, . .., ajv; h,..., bN) = |ai, . ■ . ,ak - e, . . . gn; bi,. .., 6jv) 

and 

Vk{K)\ai, . . .,aN;bi, . . . ,6jv) = |ai, . . .ajv;6i, . . . ,bk - k, ■ ■ ■ ,bN) 



(53) 



(54) 



(55) 



Thus, with the choice of Gaussian smearing functions sharply peaked at the points r/j in the operators Rj and 
the basic commutators are 



[Rk,Ui{e)] = -2lUSkiUi{^), 



(56) 



Since R^ operators have zero in the spectrum (ic. there are states such as |ai • • • , = 0, • • • a„) that have zero 
eigenvalue), there is no inverse operator and an indirect definition is required. This is achieved by Poisson bracket 
identities such as 



{^\,Ui{e)} 



first noted by Thiemann. The operator 



1 



Uki-e) 



^/S^,f>fe(e) 



(57) 



(58) 



depends on the parameter e through its dependence on the translation operators UkC. Its action on the basis states is 



— |ai,...,ajv;6i,...,6iv) = ^j^iVl^k - e| 



\ak\) |ai,...,aAr;6i,...,6jv)- 



(59) 



A symmetric version of this operator is similarly defined. 

Let us now turn to realizing operators corresponding to spatial derivatives of field variables such as (f)' and R'. To 
do this in a controlled manner it is easier to first restrict the choices of points so that they define a uniformly spaced 
lattice with spacing A. This is the first place where we explicitly introduce a radial lattice. It also has the effect of 
selecting a subspace of the full Hilbert space we defined above, by working only with a fixed set of radial points. We 
can now follow what is done in numerical methods, ie. 



(60) 



9 



and 



/ {rk) , (61) 

for any function /fe. This first of these suggests the operator 

= ^^1^' (62) 

with similar expressions for other operators. Such choices are of course not canonical, and may be viewed as an 
additional freedom in realizing a quantum theory. 

Lastly the lattice local obsevable momentum operators are 

pRk ^ ^^(Me) - Ulie)), PL = J(2 - Me) - Ulie)), (63) 

P<i>k ^ ^m^) - P^k ^ §(2 - V,{k) - V^{k)), (64) 

with action on basis states given by 

Ip 

P<t>k\a;h) = —{\a;bi,...,bk - K,...,bN) - |a; 6i, . . . , 6fe + /t, . . . , 6jv)), (65) 

and 

-P|fel«;^) = ^(2|a;6) - \a;bi, . . . ,bk - k, . . . ,bN) - |a;6i,.. .,6^ - k, . .. ,6iv)) (66) 

In constructing more complicated operators the question of the operator ordering is important. One option is a 
symmetric ordering, 

AB^AB= {Ab + bA) /2. (67) 

Another possibility is an order at which, e.g., R is to the right of Pr. We will see in the following that it has an 
advantage of annihilating the state of zero gravitational excitations. 



B. Constraint operator 

Our goal in this section to use the above definitions of basic operators to give a prescription for the reduced 
Hamiltonian (20). The main issue is the definition of the square root VY in this Hamiltonian, where 

Y = {PRRf -n- l&R^H^, n = 16R^{2RR" - 1 + R''^). (68) 

An operator representing this expression may be defined using a Dirac-like trick by suitably extending the kinematical 
Hilbert space [26]. However there is an alternative way to write the constraint (21) such that the square root problem 
is bypassed. Substituting for Pj^ from Eqn. (14) into the diffeomorphism constraint gives 

Cr = -P'nR-^+Pc^cP'^0 (69) 

If supplemented with the requirement F > (which is the case for classical solutions), this is equivalent to 

- iy' + {P^cP' - P'jiR)VY ~ 0. (70) 

This implies the constraint 

C = - {P4,(j)' - P'uRfY ~ (71) 
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In this form it does not contain a sqiiarc root and it is now straightforward to construct the corresponding operator 
using the basic ones defined above. This path for constructing the constraint operator is one of the results of this paper, 
and provides an alternative to the Dirac-like method used earher [26] where a square root operator is constructed. 

Before giving a construction of an operator analog of (71) , we note however that this constraint may have quantum 
solutions, which are not solutions of the original constraint. There is however an obvious check using Eqn. (19): after 
obtaining a solution, we must check that the expectation value of the operator analog of 



Y' - {P^cj)' - PrR') f {PrR! + P^^')dr' 
Jo 



2 



(72) 



in the proposed solution does not vanish. As we discuss below it is possible to construct the corresponding operator. 
Let us first focus on the local operator Y^. Its ingredients include 

^^|a,6)=4/|.o^|a;6), Kl\a;b) = ^ (V|afc - e| - ^^)^ \a;b), (73) 

and 

4/2 

{^if\a;b) = ^{bk+i-bkf\a;b), (74) 
and since all factors of the field Hamiltonian (density) commute, 

H^,= '^PlKl + ^Rl{4>',r. (75) 
To complete the construction of Y one needs also 

4/2 

{R',f\a;b) = j^{ak+i- akr\a;b), (76) 

and 

Rh\a;b) = -^{ak+i - 2ak + ak-i)\a;b). (77) 

The last term is 

P^R'^'^{P%R' + R'P^), (78) 

which acts as 

= -^(24|a;^)-[(afc-e)^+a^]|ai,...,afe-e, ...,ajv;6)-[(afc+e)^+4]|ai,...,afe+e, ...,ajv;6)) (79) 
Putting these pieces together we get 

Yk = (Ppi^)k - l&RlmkR'L - 1 + R'k ) - l6RlH4,k (80) 

To complete the constraint operator one also needs the commuting operators {P^(j)')k, {PRR)k- Both are obtained 
by, e. g., applying the symmetric quantization condition Eq. (67) to the elementary operators of the previous section. 

{P^(j)')k\a; b) = {bk+i -bk + ^K)\a; . . . ,bk - k, . . .) - {bk+i - bk - ^K)\a; . . . ,bk + k, . . .) j (81) 

and 

{PRR)k\a;b) = T^\^ak\...,ak+i-e,...;b) - {ttk- ^e)\...,ak-e,...;b) 

-ak\..., ak+i + e, . . . ; 6) + (ttfe + |e)| . . . , afe + e, . . . ; 6)^ (82) 
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C. Hamiltonian operator 

The asymptotic Hamiltonian is the surface term in (7) for which we need to define an operator for Pa. A direct 
quantization of the Hamiltonian density in Eqn. (20) would require us to define the operator 

^phys ^ N^'^^p^B)^ + |y,|i/2 + ((P^)k + W)k) , (83) 

which has the square root term jlfcl^/^, just as for the constraint. Since Yk is not diagonal in the basis we are using, 
this operator is not easy to define unless we go to a different basis. However, from the constraint Cr we see that 
classically on the constraint surface we have 



PA{r) = r dr'{P^cj>' + PrR'). 
Jo 



(84) 



This suggests that for physical states it is possible to compute the energy by finding an operator analog of the r.h.s. 
of this equation. Since the quantization we are using utilizes a radial lattice we can write the integral as a discrete 
sum over the lattice points r/j. It is therefore reasonable to suggest the definition 

k 

PAk = Mp J2 + (Wi)] (85) 

This operator is useful with the type of scmiclassical states we define below in Sec. V. In such states a computation of 
the expectation value (^|PAfe|V') is possible. The energy of the quantum spacetime would be this expression evaluated 
at the farthest lattice point, ie. we would take the limit fc — >■ oo after computing the expectation value. This would a 
possible analog of the classical definition where the energy 

E= lim PA{r)N''{r). (86) 

We note that the asymptotic falloff of is determined by the classical requirement of functional differentiability [8] , 
and this behavior of N'^ carries over to the quantum theory, since N'^ is not a phase space variable. 



IV. PHYSICAL STATES AND ENTANGLEMENT 



Initial states of the gravity-matter system should satisfy the quantum constraint 

C>) = 0, (87) 

which is supplemented by the requirement that \7p) belong to the positive part of the spectrum of |y|. It remains to 
be seen whether any of the operator ordering choices allows for this to be satisfied on a sufficiently large set of states, 
or if another realization of the constraint is necessary to accomplish this. Nevertheless, it is already possible to make 
a few remarks. We give here two observations, on gravity-matter entanglement and on the information loss problem. 

Firstly, in the case of pure gravity (0 = = 0), the ordering that puts the R operator to the right results in the 
operator form of the constraint (71) gives 

C|0) = (y'^/A - P'iB?Y^ |0) = 0, (88) 

where |0) stands for the state with no excitations, ie. all Ui = bi = 0. This may be viewed as a "degenerate metric 
vacuum" because it is the eigenstatc of the field operator Rk with zero eigenvalue at all points r^. 

Secondly, the presence of the matter-gravity terms in the constraint turns it into entangling operator (see eg. [23]). 
This can be seen by considering the constraint in the form of Eq. (71). It splits cleanly into gravity part, and a 
separate gravity-matter interaction 

C = Cg<»Im + J2 ® (89) 

a 

The latter term contains monomials that involve gravity and matter operators. Consider their action on the basis 
states of the kinematical Hilbert space H = 'Hg ® ^m- The term P^kPRk serves as an example: 



P^PrW) ^\b) (X {\ai- ■ ■ Qk - e- ■ ■ ajv) - |ai • • • Ofe) - e • • • ajv)) «) (|6i • • • 6fe - k • • • 6m) - \bi ■ ■ ■ bk + k- ■ ■ bm)) , (90) 
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which is a direct product of entangled gravity and matter modes. However, adding different terms in the constraint 
results in a superposition of such direct product states, and the constraint operator cannot be written as C = Cg®Cm- 
Some states of the form |V')g <H) |<^)m may be transformed into a direct product form, 

C{\i^)G\^)M) = W)GW)M, (91) 

provided that they satisfy an additional constraint 

^C^^Cf (|^)g|¥')m)(x |V')Gk)M. (92) 

a 

As a result, even if such states exist, they form a lower-dimensional set then the set of all physical states, and 
by themselves do not form a vector space: their linear combinations are by definition entangled. Since we expect 
the semiclassical states to be suitably defined coherent states, we see that these configurations are geometry-matter 
entangled. This result of course holds regardless of weather a state corresponds to a classical black hole. The 
entanglement prevents the direct product gravity— matter decomposition of a generic physical state. 



V. SEMICLASSICAL STATES 



Semiclassical states that are peaked at a given classical configiirations may be defined for a field theory just as for 
a quantum mechanical systems. For a constrained system such as the one we are considering, it is possible to obtain 
states that are peaked on classical constraint free data. These may be viewed as approximate physical states in a 
precise sense. 

For a constraint C in a field theory, the steps we propose are as follows: (i) begin with an explicit classical solution 
of the constraints, such as that given in Sec. II. B, (ii) select sampling points Tj (or a lattice) , and the corresponding 
solution values R{ri), 4'{ri) etc. , (iii) construct a Gaussian state peaked at each field value i?(r,) and (jiiri), and 
finally (iv) take the product of the Gaussians, one at each point to give the field theory semiclassical state. This is 
what we describe in detail below. 

Reasonable requirements for an approximate semiclassical state are that 

{C{rk) = 0, {C{rkf) = 0, (93) 

for each sample point r^. We give here a generalization of the semiclassical states for Friedmann- Robertson- Walker 
cosmology given in [? ], and show how this can be utilized for the present model. 

Let us consider first a single lattice point and the basis states at this point defined by \ak) = \me)k, where m in 
an integer. Let us consider the state at rk defined by the linear combination 

oo 

\P°h,RX': = ^ E e-t(-)^e™^«°e^-^^«|me),,. (94) 

m=— oo 

This is a Gaussian state of width t (measured in Planck units), where the (real) parameters P]2. and i?° are the field 
values corresponding to a classical configuration at the point rk- The width t is a measure of how strongly peaked 
the state is on a given classical configuration, ie. t <^ 1 means that fluctuations around P^,BP are small. We shall 
see in the following that this state has a number of desirable properties, and because of this, it allows a construction 
of approximate physical states of the theory we are considering. (There may be other non-Gaussian states with such 
properties but our purpose here is to point out that this one is useful.) 

The normalization condition for the state gives an expression for the constant C. It is the convergent sum 

e-*^'™'e2«°^™. (95) 

m=— oo 

Calculation with this state gives the expectation value [? ] 

{U^) = e'^-PR('-'=) e-*^'/^ K{e, t, R°). (96) 

where 

,' l + 2E^^oCos[^(l + ^)]e-^-V*' ^ 
K{e, t,R)-\ ^ ^ ^ ^^^^ ,_..„Vte^ • (^^^ 
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Equation (96) together with the definition (48) gives the expectation value 

^ ^ sin[pO^(.,)e] i^(e,t,J?°). (98) 

This formula has the limits 

lim (P^) = sin(pO(rfe)e)/e, (99) 
lim (P^) = pO(rfc). (100) 

The first shows that the semiclassical state in field space is peaked at the corresponding phase space value. The 
second shows that the field continuum limit of the momentum expectation value has the appropriate peaked value in 
this state, even though only field translation operators exist in the representation we are using. 

These semiclassical states defined at each point r, can now be used to give a state for the entire radial lattice given 
any classical field configuration R{r),Pii{r) by taking a product of the point states over the lattice {rk}k=n ie. 

N 

\R{r), P{r)Y'^ := [] \R{rk), PR{ru))Y, (101) 
fe=i 

The simplest semiclassical state for both sets of fields is the product 

\x)^\R{r),P{r)Y^^ W),P^{r)Y^'. (102) 

It is not difficult to configure entangled products of semiclassical states, given two distinct classical solutions of the 
constraints. 

Given that the expectation values in these states give the classical peaking values it follows that the states peaked 
on classical constraint free data, such as those constructed in Section II, satisfy the expectation value 

{x\C\x)=^ + 0{n, (103) 

where cr > 1. We note also that the constriiction we have given is quite specific, and there is control on the fluctuation 
(C^), which is itself a function of the width t of the point Gaussian states. Its value can be tuned to reduce the 
fluctuations. It is in this sense that these states are approximately physical for sharply peaked states t One can 
now compute quantities such as the expectation value of the energy in such states, which is tedious put possible. 

A. Effective equations 

One approach for obtaining quantum gravity corrections is to use semiclassical states to derive "effective con- 
straints." These can then be used to obtain modified equations of motion in the usual way. This approach is 
qualitatively related to the ideas underlying Ehrenfest's theorem in quantum mechanics. One computes expecta- 
tion values of the constraints in states peaked on classical configurations that are not solutions of the initial value 
constraints to obtain 

(<^,P^,i?,PK|(7|</.,P^,i?,PK) = C^''^^^^='^'(</.,P^,i?,Pfl)+0(r). (104) 

It is understood here that the classical term on the right hand side is a function of the phase space configuration on 
which the states are peaked. The corrections are functions of the state's width t in Planck units. 

This approach has been used for homogeneous isotropic cosmology [28], but requires careful scrutiny in the field 
theory case here. One has to check that the Poisson algebra of the effective constraints closes with the state width 
corrections included, since these corrections are functions of the phase space peaking values of the semiclassical state. 
One of the goals of the program presented here is to obtain consistent effective equations of motion starting from 
semiclassical states. 



VI. SUMMARY AND DISCUSSION 

The paper contains a number of developments beyond earlier works on the spherically symmetric gravity-scalar 
model: (i) A construction of semiclassical states for the quantum theory of gravity coupled to a minimally coupled 
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masslcss scalar field. These states may be utilized for constructing scmiclassical effective constraints, which may then 
be used for classical numerical evolution. If used in conjunction with constraint free data, such as that given in Sec. 
II. B, this construction gives the first known examples of physical semiclassical states, (ii) An alternative construction 
of the Hamiltonian operator in a fixed time gauge which can serve as a starting point for Monte Carlo simulations of 
the quantum theory. Since this model is a two-dimensional system that is effectively written as a lattice theory, this 
method has the potential to reveal interesting non-perturbative quantum phenomena such as phase transitions. This 
is at present being studied. 

The form of the Hamiltonian also reveals that matter-gravity entanglement is an inherent feature of evolution in 
quantum gravity; the action of the Hamiltonian on a product state gives an entangled state after a single evolution 
step. In the full quantum problem the constraint operator forces a bipartite entanglement (with the two subsystems 
being the kinematical Hilbert spaces of matter and gravity degrees of freedom, respectively) already on the initial 
physical states. It is tempting to summarize this observation by modifying Wheeler's one line description of general 
relativity [29]: quantum gravity tells geometry and matter how to entangle. 

The degree of entanglement may be computed in the usual way by tracing a density matrix over cither matter 
of geometry degrees of freedom. If a state describes a black hole there is a second entropy that may be calculated, 
namely that obtained by tracing the (pure) density matrix over both geometry and matter degrees of freedom inside a 
trapped region. This of course would be different from the usual entanglement entropy where the trace in the interior 
applies only to the matter degrees of freedom. 

It will be apparent to the reader that the extraction of physical results from this formalism requires further work. 
This could proceed along two distinct avenues. The first is the derivation of scmiclassical effective equations where 
one would compute expectation values of the constraint operators in suitable states to obtain "effective" constraints. 
These would then be used as the basis of a quantum gravity corrected classical dynamics which could be integrated 
using numerical methods. Among other things it would be useful to sec what becomes of the critical scaling observed 
at the onset of black hole formation, especially the critical solution. Initial results using this approach [30, 31] indicate 
that black holes form with a gap, but so far there are no results on the critical solution itself; with expected singularity 
avoidance, it is likely that this is replaced by a critical and unstable bound state - a finely tuned boson star. 

The second approach to studying dynamics could utilize Monte Carlo methods. The quantum model as formulated 
here is like a statistical mechanical system, with the difference that it is still a constrained theory. One can imagine 
sampling from phase space such that the samples are solutions of the constraints up to some threshold, along with the 
usual Monte Carlo selection criteria. Such an approach would have the potential to yield fundamental information 
concerning phase transitions. 
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